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SUMMARY 


PROBLEM 

Detection  of  a sinusoid  obscured  by  white  noise. 

RESULTS 

1.  Currently  employed  methods  for  the  detection  of  a sinusoid  obscured  by  white 
noise  were  related  to  three  classical  problems  of  detection  theory.  The  discrete  Fourier 
transform,  the  power  spectrum,  and  the  averaged  power  spectrum  each  were  shown  to  be 
the  optimal  processor  (iti  the  sense  of  making  a least-risk  decision)  for  a particular  signal 
model. 

2.  The  first-  and  second-order  statistics  under  H|  (signal  plus  noise)  and  Hq  (noise 
alone)  were  compiled  for  seve  al  related  processing  structures  of  current  interest.  Included 
were  the  autocorrelation  function  and  the  recursive  exponential  correlator. 

RECOMMENDATIONS 

1.  Extend  the  ROC  (receiver  operating  characteristic)  performance  results  of  this 
report  to  include  the  autocorrelation  function  and  recursive  exponential  correla:  or  process- 
ing structures. 

2.  Investigate  performance  via  the  ROC  curve  of  the  discrete  Fourier  transform, 
autocorrelation  function,  and  recursive  exponential  correlator  processing  structures  when 
the  sinusoid’s  frequency  varies  over  the  observation  interval. 
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PROCESSOR  STRUCTURES  FOR  THE  DETECTION  OF  A 
SINUSOID  OBSCURED  BY  WHITE  NOISE 

I.  INTRODUCTION 

This  report  has  two  main  thrusts.  Piimarily,  the  desire  is  to  relate  currently 
employed  methods  for  the  detection  of  a sinusoid  obscured  by  white  noise  to  three  classical 
problems  of  detection  theory.  It  will  be  shown  that  the  discrete  Fourier  transform,  the 
power  spectrum,  and  the  averaged  power  spectrum  are  each  the  optimal  processor  (in  the 
sense  of  making  a least-risk  decision)  for  a particular  signal  model.  The  utilization  of  one  of 
these  processors  implicity  assumes  that  particular  signal  model  to  exist  and  thus  defines  an 
upper  bound  on  performance  even  if  the  actual  received  signal  process  could  admit  a proc- 
essing scheme  whose  performance  would  be  superior.  Of  secondary  importance,  the  first  - 
and  second-order  statistics  under  H|  (signal  plus  noise)  and  Hq  (noise  alone)  are  compiled 
for  several  related  processing  schemes  of  current  interest.  Included  are  the  autocorrelation 
function  and  the  recursive  exponential  correlator.  These  statistics  will  prove  useful  in  pre- 
dicting the  performance  of  such  schemes  in  light  of  the  processor  structures  already 
mentioned. 

Our  fundamental  goal  is  the  detection  of  a sinusoid  obscured  by  white  noise.  If  the 
form  of  the  received  signal  was  known  exactly,  the  optimal  detector  would  simply  threshold 
the  output  of  a filter  matched  to  the  signal  waveform.  Unfortunately,  the  actual  situation 
encountered  is  one  where  phase,  amplitude,  and  frequency  uncertainty  usually  exists. 

A commonly  employed  processor  structure  for  detection  is  a display  of  the  discrete 
Fourier  transform  magnitude  squared  (|DFT|“)of  the  observed  time  series  (i.e..  the  power 
spectrum).  In  addition,  averaging  of  the  power  spectra  from  several  consecutive  data  sets 
often  is  performed  either  visually  or  automatically  to  increase  the  detectability  of  particu- 
larly weak  signals. 

This  report  will  discuss  several  potential  configurations  for  the  detector  structure. 
The  reference  point  for  their  comparison  will  be  the  following  two  hypothesis  testing 
problem: 

H| : X(n)  - S(n)  + n(n) 

Hq:  X(n)  = n(n).  ( l.’J 

where: 

0<  n < N-l 

x(n)  the  observed  sequence 

S(n)  = A cos  ( 2irf0nA  + 0)  the  signal  sequence 

n(n)  the  white  Gaussian  noise  sequence;  n(n)  ^ N(0.o“). 


| 
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The  assumption  is  made  that  the  observed  sequence  was  obtained  from  a continuous  time 
series  (0  < t < (N-l)A)  bandlimited  to  W hertz  and  sampled  every  A = 1/2W  seconds. 

The  motivation  for  considering  several  detector  structures  can  be  found  in  the  details 
of  any  one  particular  problem.  Due  to  considerations  of  computational  speed  or  the  under- 
lying physics  of  a given  situation,  one  structure  may  be  preferred  over  another.  For  example, 
discrete  Fourier  transforms  are  faster  to  compute  than  autocorrelation  functions.  However, 
the  ability  to  window  the  autocorrelation  function  before  transforming  to  the  frequency 
domain  provides  a certain  degree  of  flexibility  in  what  the  resulting  display  looks  like. 
Another  concern  is  stability  of  the  sinusoid.  If  its  phase  or  frequency  vary  with  time, 
averaging  of  power  spectra  computed  from  short  consecutive  data  sets  may  be  necessary  to 
achieve  the  best  detectability. 

The  organization  of  this  report  is  as  follows.  In  Section  II,  the  likelihood  ratio  proc- 
essor is  defined  and  its  performance  evaluation  in  terms  of  the  ROC  (receiver  operating 
characteristic)  curve  is  discussed.  The  next  three  sections  develop  three  optimal  receiver 
structures  - SKE  (signal  known  exactly),  SKEP  (signal  known  except  for  phase),  and  SKEP 
Independent  Increments.  Section  VI  defines  the  discrete  Fourier  Transform  (DFT)  and 
relates  it  to  each  of  the  three  optimal  detector  structures  mentioned  above.  The  autocorre- 
lation function  and  its  relationship  to  ti  e power  spectrum  of  a sequence  are  the  subject  of 
the  following  section.  Lastly,  Section  VIII  develops  expressions  for  the  REC  (recursive 
exponential  correlator). 

II.  DETECTION  THEORY 

A.  The  Likelihood  Ratio 


A helpful  baseline  from  which  to  compute  various  detector  structures  is  that  of  the 
likelihood  ratio  processor.  Let  the  vector  of  observables  be  denoted  by 

X = (X(0) X(N-1)IT.  (2.1) 

Based  upon  the  observation  vector  X e x.  the  processor  must  make  a decision  (Dq  or  D] ) as 
to  which  hypothesis  it  believes  is  true.  Classical  detection  theory  has  shown  that  decisions 
based  upon  the  likelihood  ratio  are  optimum  for  a wide  range  of  goodness  criteria  (Peterson. 
Birdsall,  and  Fox,  1954) 


A(X)  d 


p(X|H,) 

P(X|H0) 


D| 

^0 


r\- 


(2.2) 


Birdsall  has  shown  more  generally  that  any  optimality  criterion  based  on  “detection 
probability”  P(D  j |H  j ) and  “false  alarm  probability”  P(D|  |Ho)  where  good  decisions  are 
preferred  over  bad  leads  to  the  calculation  of  A(X)  as  the  decision  statistic  (Birdsall.  1973). 
Thus,  a separation  is  achieved  between  the  processing  of  X and  the  actual  optimality  crite- 
rion chosen  which  arises  in  the  selection  of  a threshold  value  r?. 

The  situation  may  arise  where  one  or  several  parameters  under  either  or  both  Hq  and 
H|  are  uncertain.  These  are  then  modeled  as  random  variables  and  any  prior  knowledge 
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about  them  is  summarized  by  a priori  probability  density  functions  p(0q)  and  p(0  | ).  The 
desired  decision  statistic  now  becomes  the  ratio  of  marginal  probability  density  functions 
on  x: 


A(X) 


P(X|0,,H,)p(0,)d0, 


P(X|£q,  Hq)  p(  £ q)  ^-0 


(2.3) 


where  £j  e0|  and  Oq  cOp.  A block  diagram  of  the  processor  structure  implied  by  (2.3)  is 
given  in  Figure  2.1. 

B.  Performance 


The  complete  description  of  a detection  device  includes  both  the  processor  itself 
(i.e.,  the  mathematical  transformation  from  observation  space  to  decision  statistic)  and  the 
performance  of  the  processor  evaluated  with  respect  to  the  goodness  criterion  initially 
chosen.  As  mentioned  earlier,  the  likelihood  ratio  has  been  shown  optimum  for  any  good- 
ness criterion  based  on  “detection  probability”  P(D  | |H  j)  and  “false  alarm  probability 
P(D  i |Hq)  where  good  decisions  are  preferred  over  bad.  Thus,  the  appropriate  description 
of  performance  for  a likelihood  ratio  computing  device  is  its  detection  and  false  alarm 
probabilities  as  a function  of  decision  threshold.  The  precise  definition  of  these  terms 
(which  arise  from  within  a RADAR  and  SONAR  context)  now  will  be  given. 

Since  the  likelihood  ratio  is  simply  a transformation  of  random  variables  (the  obser- 
vation vector  X)  to  a one-dimensional  decision  statistic  (A(X»,  the  likelihood  ratio  itself 
will  be  a random  variable  whose  probability  density  function  will  depend  on  which  hypoth- 
esis (Hq  or  H | ) is  actually  active  on  x-  Recalling  that  the  threshold  rj  divides  the  decision 
space,  define 

PD  £ PtDjIH,)  d / p(A|H,)dA  (2.4) 


\< 


Pto,)  P<V 


Figure  2.1.  Likelihood  ratio  processor. 
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PF  £ P(D,|H0)  * 


*/ 


p(A|H0)  dA. 


(2.5) 


Peterson,  Birdsall,  and  Fox  introduced  a graphical  representation  of  Pj)  versus  Pp  as  a func- 
tion of  tj  known  as  the  ROC  (receiver  operating  characteristic)  curve  (Peterson,  Birdsall,  and 
Fox,  1954).  The  ROC  curve  will  be  the  means  by  which  performance  of  the  detection 
receivers  discussed  in  this  report  will  be  evaluated  and  compared. 

III.  SKE  (Signal  Known  Exactly) 

In  the  SKE  problem,  the  received  signal  (when  present)  is  known  exactly.  The 
statistics  are  Gaussian  under  both  H j and  Hq. 


A(X)  = 


p(X|H ! ) 
P(X|H0) 


(the  likelihood  ratio) 


1 

■ N-i 

- £ (X(n)-S(n))2 
n=0 

(2tto2)N/2  P 

. 2o2 

1 

■ N-l 

- £ x<n>2 
n=0 

(2»o2)N/2  CXP 

. 2o2 

= exp 


= exp 


N-l 


X(n)  S(N)  S(rA- 


2 o- 


n=0  °~ 

/N-l  N-l 

£ X(n)S(n)-|^  S(n)“j 
” n=0 


r.=0 


Note:  Bandlimited  white  Gaussian  Noise 

^0/~ 


a-  = N0W 


W Hz 


(3.1) 


E ^ signal  energy 


(3.2) 
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Thus, 


n=0 


N-l 

O’Vo 


A(X)  = exp 


N0 

N-l 


-L  S X(n)S(n)--|- 
n2  Na 


1°"  n“0 

Any  monotonic  function  of  A(X)  will  also  be  an  optimal  test  statistic. 

N-l 

CnA(X)  = -~  + 4;  7 X(n)  S(n) 

N0  °2  n=0 


/•NA 

-*+±[  5 

N0  NoJo 


X(t)  S(t)  dt. 


(3.3) 

(3.4) 


(3.5) 


(3.6) 


The  last  term  in  each  of  the  above  two  expressions  is  known  as  a “matched  filter.” 
Substituting  the  expression  for  S(n), 


N-l 

^ X(n)  cos(27rfQnA) 
n=0 


■ a( 


En  A(X)  + 


(3.7) 


Since  the  parameters  A,  o~,  E,  and  Nq  are  known  in  the  SKE  problem,  the  above  is  simply  a 
monotonic  function  of  A(X).  Thus,  a = Re{X(k)},  k = NfgA  is  an  optimal  decision  statistic 
|see  Section  VI  for  the  definition  of  X(k)l . Computing  the  DFT  of  the  observed  sequence 
wiil  be  the  optimal  processor  when  the  signal  is  known  exactly. 


E|a|H,  | 

= 4n 

E|a|H0l 

= 0 

var|a|Hj  1 

= var|a|H0l  = 77 

In  general,  the  entire  ROC  curve  is  necessary  to  completely  specify  performance. 
However,  in  the  SKE  problem,  performance  is  summarized  by  a single  number  known  as  the 
detectability  index  d-.  In  this  case,  the  distribution  of  the  optimal  decision  statistic  is 
Gaussian  under  H | and  Hq  with  equal  variances.  By  definition  (Van  Trees,  1968) 
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2 A (E(a|H!  ] -ElalHpi)- 


d-  6 


var|a|H0l 


A“N 
2 o“ 


2E 

N, 


|=N*SNR;SNR  = A_ 


2 o“  ^ 


(3.8) 


0 


Detection  and  false  alarm  probability  expressions  corresponding  to  Eqs.  (2.4)  and  (2.S)  are 

(3.9) 


-.(^4 


Pp  = erfc 


(3.10) 


where 


*00 

erfc-  ' / ^jexp(f) 

•'x 


dx. 


(3.11) 


The  SKE  performance  curves  are  illustrated  in  Figure  3.1  on  normal-normal  paper.  Note 
that  performance  increases  linearly  on  the  negative  diagonal  as  a function  of  d. 

IV.  SKEP  (Signal  Known  Except  for  Phase) 

In  the  SKEP  problem,  the  received  signal  (when  present)  is  known  except  for  its 

phase: 

S(n)  = A cos  (2jrfgA  + 0), 

where:  0 is  a random  variable  distributed  uniformly  between  -jt  and  it. 

The  statistics  are  conditionally  Gaussian  under  H | and  Gaussian  under  Hq. 


■/' 

■if 

J-Tt 


A(X|0)  p(0)d0 


exp 


N-l 

- T7-  + A T X(n)  cos  ( 2fffnnA  + 0) 

No  ,(=0 


d0 


= exp 


= exp 


U1  T-  f A<“2  + cos  («  + 01 00 

L NoJ  -*L,  0- 


(4.1) 
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where: 


X(k)  = a + ib 
k = Nf0A 

d = tan-'  (b/a). 

Note  the  following  monotonic  function  of  the  likelihood  ratio: 

A(X)  « (a“  + b“)  = |X(k)|",  k = Nf0A.  (4.2) 

Thus,  a-  + b-  = |X(k)|~  is  an  optimal  decision  statistic  |see  Section  VI  for  definition  of 
X(k)J . Computing  the  power  spectrum  of  the  observed  sequence  is  the  optimal  processor 
when  the  signal  is  known  except  for  phase,  with  the  phase  having  a uniform  distribution. 

The  ROC  curves  of  the  optimal  detector  for  the  SKEP  problem  are  given  in  Figure 
4.1  (Roberts,  1965).  The  loss  in  performance  due  to  phase  uncertainty  can  be  seen  by  com- 
paring the  SKEP  curves  with  the  SKE  curves  for  the  same  value  of  2E/Nq.  (See  also  Van 
Trees,  1968  ) 


Two  important  points  must  be  stressed  with  respect  to  this  problem.  First,  it  appears 
that  the  loss  in  performance  from  SKE  can  be  summarized  by  noting  that  the  test  statistic 
contains  twice  as  much  noise  (i.e.,  simply  a 3-dB  loss).  Using  this  reasoning,  the  SKEP  ROC 
curve  in  Figure  4. 1 for  2E/Ng  = 9 roughly  should  be  parallel  to  the  SKE  ROC  curve  for 
2E/Nq  = 4.  Note  that  the  two  ROC’s  cross  each  other,  indicating  a more  fundamental  dif- 
ferential in  performance  than  can  be  accounted  for  by  a simple  adjustment  of  2E/Nq. 


The  second  point  concerns  the  use  of  a detectability  index  to  compare  performance 
between  two  processors.  A detectability  index  for  the  SKEP  processor  can  be  defined  as  in 
Eq.  (3.8). 


I^(k)  = ^j-|X(k)|“,  k = NfgA  (See  Section  VI.) 

= I(a2  + b2) 

N 

E|IN(k)|H0]  = o: 


E(IN(k)|H,|  = o“  + — N 


- „4 


var  UN(k)|H0l  = o 


var  (lN(k)|H|  1 = 04  + o~^-N. 
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Define: 


'f 


2 A (E[ lN(k)|H 1 1 -E[lN(k)|H0))2 
dSKEP  “ var  [iN(k)|H0l 


On  the  basis  of  Eqs.  (3.8)  and  (4.3)  alone,  one  might  incorrectly  conclude  that  the  SKEP 
processor  will  perform  much  better  than  the  SKE  processor  when  E/Nq  > 2.  Recall  that  the 
only  valid  performance  comparison  is  in  terms  of  the  ROC  curve  as  shown  in  Figure  4.1. 

(In  retrospect,  it  can  be  seen  that  Eq.  (4.3)  does  not  take  into  account  that 

var  [l]iq(k)|H  ] ] var  [I]s](k)iHol  and  that  the  distributions  of  the  test  statistic  are  now 

chi-square  instead  of  Gaussian). 


V.  SKEP  (Signal  Known  Except  for  Phase)  INDEPENDENT  INCREMENTS 

In  the  SKEP  Independent  Increments  problem  the  received  signal  (when  present)  is 
known  except  for  its  phase 

S(n)  = A cos(2rrfQnA  + 0g), 


where: 


0g  is  a random  variable  distributed  uniformly  between  -v  and  ir. 


0 j is  constant  for  ?*N  < n < (C+1)N,  8 = 0,...,  L-l 

but  changes  independently  between  increments. 


Data  sequence  consists  of  L blocks,  each  of  length  N.  The  statistics  are  condition 
ally  Gaussian  under  Hj  and  Gaussian  under  Hq. 


where: 


L-l 

A(X)  = PI  A(Xj)) 
8=0 


(5.1) 


Xc  = [X(fc-N) X((K+1) -N- 1)1T 

X^(k)  = ;i{)  + ibg.  k = NIqA. 
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Note  the  following  monotonic  function  of  the  likelihood  ratio: 


' 8nA,X)+f  • 


Furthermore, 


fin  l0(x)  3s  x - 8nv2irx, 
*> 

“ 4 ’ 


x>>  1 

X < < 1. 


When 


N*SNR  < < 1 
L-l 


h ■ ($] 


6-0 

Thus,  averaging  power  spectra 
L-l 


s?0nM 


(5.2) 


(5.3) 


(5.4) 


approximately  yields  an  optimal  decision  statistic  for  the  SKEP  independent  increments 
problem  when  N*SNR  < < I (See  Section  VI). 


Performance  of  the  SKEP  independent  increments  processor  is  given  in  Figure  5.1 
(Thompson,  1972).  ROC  curves  for  I,  5,  10,  and  20  independent  increments  are  illustrated. 
The  SKE  ROC  curve  for  the  same  2E/Nq  value  is  included  as  a reference.  Note  that  L=l 
corresponds  to  the  SKEP  processor  discussed  in  Section  IV.  It  can  be  seen  easily  from 
Figure  5.1  the  dramatic  effect  signal  instability  has  on  detectability.  (Performance  for  the 
detector  structure  in  Eq.  (5.4)  for  arbitrary  SNR  can  be  found  in  Marcum,  1960.1 

VI.  DISCRETE  FOURIER  TRANSFORM 

A.  Discrete  Fourier  Transform  (DFT) 


N-l 

X(k)  A ^ X(n)  e"<:^N)nk 
n=0 

= a + ib;  a,b  ~ N |o,-^  o“j  under  Hq. 

X(k)X(k)*  = a-  + b- 
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E(X(k)|H,)  = 4N:  k = NfgA 

" 

= 0;  k # NfgA 

var  [X(k)|H0]  = No2 

4 

, i 

var  [X(k)|H,  1 = No2 

Computing  the  DFT  of  the  observed  sequence  is  the  optimal  processor  when  the  signal  is 
known  exactly  (See  Section  III). 

B.  Power  Spectrum  (Periodogram) 

IN(k)  $ -jJj  |X(k)|2 

= ^ (a2  + b2);  a,  b,  ~ N |o, -j o2j  under  H0. 

(See  also  Oppenheim  and  Schafer,  1975.) 

E(IN(k)|H0l  = o2 

i A2 

E[IN(k)|H|l  = o2+^-N;  k = NfgA 

= o2:  k * NfgA 

var  (lN(k)|Hgl  = o4 

var  (lN(k)|Hj ] =a4  + o2~N;  k = NfgA 
= o4:  k * NfgA. 

Note  that  the  statistics  of  Ij^(k)  are  chi-square.  Computing  the  power  spectrum  of 
the  observed  sequence  is  the  optimal  processor  when  the  signal  is  known  except  for  phase 
with  the  phase  having  a uniform  distribution  (See  Section  IV). 
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C.  Avenge  Power  Spectrum 

Average  K spectra  each  derived  from  a time  sequence  of  length  N. 

K-l 

ItyOOgyg  * v ^ ^N^i 

i-0 

i-0 

K-l 

. i y i (tp + b?> 

1-0 

E[IN(k)av|lH0]  » o2 

E[IN(k)tvg|H,  1 - o2  + ^ N;  k » Nf^ 

- a2\  k # Nf0A 

var  (IN(k)avg|H0)  - 

v«r(lN(k)IV,IH,|  - i |o4+c2^n};  k - NfoA 
“ £o4;  k # NfgA. 

Avenging  power  spectra  computed  from  consecutive  blocks  of  the  observed  sequence  closely 
approximates  the  optimal  processor  for  N-SNR  < < 1 when  the  signal  is  known  except  for 
phase,  with  the  phase  being  uniformly  distributed  in  each  increment  and  independent  between 

increments  (See  Section  V). 

r 

VII.  AUTOCORRELATION  FUNCTION 
A.  Autocorrelation  Function 
N-l-m 

Rxx(m)  - ■jjj  ^ X(n)  X(n+m),  m=0, 1 , N-l . 

n*0 

Note:  Defining  Rxx(-m)  = Rxx(m) 
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«*(*.*■  *7*  ****»«*i<I»  IfrvrWlHI^  ■SMfj  }*^t»MV'l»i^r-(«'  **",<0  ■‘WT'VWMHr*^ 


IN(k)  ■ ,j“IX(k)|“ 

N-l 

= ][  Rxx(m)  e-j(2ir/N>mk 

m=-(N-l) 

(See  Section  V!  for  definition  of  1^00.) 

E(Rxx(m)|Ho]  * a 2,  m=0 

= 0,  m >0 


, A2 

E(Rxx(m)|Hj)  = o2  + ^-,  m=0 

N-m  A2  c k 
= — — cos  27rf0mA, 

var  IRxx(m)|H0l  = ^-o4,  m=0 

N-m  ^ „ 

= — — . m>0 

N2 


rn  >0 


var  [Rxx(m)|H|l  = ? 

Note  that  the  statistics  of  Rxx(m)  are  not  Gaussian. 
Define: 


Sxx(k’  5)  ~ 


L-l 

£ Rxx(m  + 5)e-j(2jr/L>mk, 
m=0 


I <L<N-6- 


Sxx(k,  5)  is  the  DFT  of  an  L-point  segment  of  the  autocorrelation  function. 


E(Sxx(k,  5)|Hq] 

■» 

= 0“, 

6 = 0 

= 0. 

6 > 0 

E|Sxx<k,5)|H|| 

■> 

= o~  + 

A2  (l  L2  \ 
2 \2  ~ 4N/ 

6 = 0:k  = Nf0A 

<h 

it 

L L:  + 2L6\ 

2 4N  /' 

5 > 0:  k = Nf0A 
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var  [Sxx(k,  5)|H01  = j^.  + J_  (-L2  + 2NL  + L)J  , 6-0 

= o4l-  - (-L2  + 2L(N  - 6)  + L)1  , 6 > 0 

l2M  J 

var  (Sxx(k,  5)|H,  ] = ? 

In  general,  Sxx(k,  6)  is  complex.  When  6 = 0,  the  expression  for  variance  yields 
twice  the  value  obtained  in  Section  VI  B.  In  practice,  only  Re|Sxx(k,  0)|  would  be  used. 

B.  Averaged  Autocorrelation  Function 

Average  K autocorrelation  functions,  each  derived  from  a time  sequence  of  length  N. 


K-l 

Rxx(m)avg  “ £ ^ ^xx^i 
i=0 


K-l  N-l-m 

= K S N £ X(n>i  X(n+m)i’  m=0’  • • • ’ N'1 

i=0  n=0 

E(Sxx(k,  6)|H0]  = o2,  5 = 0 


= 0, 


6 >0 


5 = 0;  k = Nf0A 

6 > 0;  k = Nf0A 


. AZ(  l L2  + JL5\ 

2 \2 ' 4N  / ’ 

var  [Sxx(k.  8)IH„  * £ ft ' + (-L-  + + t)]  , 

= T|^(-L:  + :L,N-s,tL)]' 


6 = 0 


6 > 0 


s 9 


var  (Sxx(k,6)|Hi 

C.  Windowed  Autocorrelation  Function 


Consider  the  autocorrelation  function  created  from  a time  sequence  of  length  N. 
Using  a triangular  sequence  (Bartlett  window),  window  this  function  back  to  length  L + 5. 
That  is. 
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Define: 


Rx?(m)  = 


• Rxx(m) 


The  above  is  equivalent  to  defining 

N-l-m 

Rx^(m)  ^ k±i_liD  i £ X(n)  X(n  + m),  m = 0, 1 , . . . , L 


E[Rx^(m)|H0)  = o^, 


n=0 

= m = 0 


= 0,  m > 0 


E[Rx^(m)|H,l  = o2+Al,  m = 0 


_ + 6 - m\  /N  - m\  A2  <,  t A 

-(-TTr-)hr ]Tcos2,,tomA- 

= 1 d 


var  [Rx^(m)|H0]  = ^ o4,  m = 0 


„ /L  + S - m\  2 /N  - m\  4 
\ L + 6 J \ n2  ) ’ 


m > 0 


var[R“(m)|H,]  = ? 


L-l 


sx^(k.  6)  a \ Rx^(m  + 6)  e"-i(2w/L)mk 

m=0 

ElS^tk,  6)|H0)  = a2, 

0 

n 

<0 

= 0, 

5 >0 

E[Sx^(k,6)|H,l  = o2+Al 

(*-«■ 

6 = 0:  k = Nf0A 

<1 

II 

l2  + :sl\ 

5 > 0;  k = Nf0A 

2 \: 

4 (L  8)/  ’ 

*6-1. 


m > 0 
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*ar  |S  “<k.  8)IH0|  = „4ii  |2 +|  + J.J  , 

L 

.o'!— ti_i|:  + 3t  '|. 

(L+6):N6  * L L-f 


8 = 0;N>>L  + 8 
6 > 0; 

N>>  L + 8 


var  (S“(k,8)|H,]  = ? 

Note  that  windowing  (VII  0)  and  averaging  (VII  B)  are  not  equivalent.  Averaging 
autocorrelation  functions  implies  that  the  data  has  been  broken  into  K independent  blocks, 
each  of  length  N.  A windowed  autocorrelation  function  is  a moving  average  function  across 
the  entire  data  sequence.  Thus,  each  data  sample  is  allowed  to  interact  with  data  up  to 
L + 8 lags  into  the  past.  In  averaging,  interaction  between  data  samples  is  restricted  to 
within  the  individual  data  blocks. 

VIII.  REC  (Recursive  Exponential  Correlator) 

Weight  vector 

W(j)  £ [W0(j),...,WL_1(j)lT 

Wm0  + I)  = (I  -0)Wm(j)  + 0X(j-8-m)X(j) 


J-l 


Wm(j)  = (I  -0)H“k  0X(k  - 8 - m)  X(k): 

k=8+m 


Wm(0)  = 0 


J[WmO)IH0]  = o2  [1 -(1 , 6 = m = 0 

= 0,  5 = 0,  m > 0 

= 0,  8 >0 

ElWm(j)IHi)  = —cos  2jrf0(8  + m)  [l -(I -0)M*+»nH]t  5>0 


var  lWm(j)|H0l 


4 *2  P - <l  - g):(H5+mHy 


_ 4/3**  1 

= cr  p-  - 


20-0- 


■1 


8 > 0 


var|WmU)|H,]  = ? 


Define: 


L-l 

Hj(k)  £ £ WmU)e^2'/L),nk 
m=0 

f 

[ 


20 


6 = 0 


E(Hj(k)|H0]  = o2|l-(l-flHl, 

= 0,  6 > 0 


E(Hj(k)|H, 


o ' jl-d-#1}  + y (t"C)  • 


6=0 


6 > 0 


where: 


C 3:  ( 1 -0>i 


-1-6  1 

2fnA 


f 1 - ( 1 - /3)~L 

Ll  -0  -0r,/foA. 


o4fl2 

var  (H:(k)|H0j  = 

J 20-02 

var  [Hj(k)|H| ) = ? 


■ l -,1-0)20-1-0  /i-n-dr211--1*' 
\ 1 -(1  -0)‘2 


» 


The  REC  is  similar  to  the  windowed  autocorrelation  function  of  Section  VII  C in 
that  both  are  moving  average  filters.  In  the  former,  the  past  is  exponentially  weighted 
instead  of  uniformly  as  in  the  latter.  Additionally,  the  REC  closely  resembles  the  weight 
vector  of  the  ALE  (adaptive  line  enhancer)  (Widrow,  et  al.,  1975)  under  the  condition  of 
low  signal-to-noise  ratio. 


IX.  SUMMARY 


The  primary  emphasis  in  this  report  has  been  the  relationship  between  currently 
employed  methods  for  the  detection  of  a sinusoid  obscured  by  white  noise  and  three  classi- 
cal problems  of  detection  theory.  It  was  shown  that  the  discrete  Fourier  transform,  the 
power  spectrum,  and  the  averaged  power  spectrum  are  each  the  optimal  processor  (in  the 
sense  of  making  a least-risk  decision)  for  a particular  signal  model.  Of  secondary  importance, 
the  first-  and  second-order  statistics  under  H|  (signal  plus  noise)  and  Hq  (noise  alone)  were 
compiled  for  several  related  processing  schemes  of  current  interest. 
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